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Abstract —The Compute-and-Forward relaying strategy 
achieves high computation rates by decoding linear combinations 
of transmitted messages at intermediate relays. However, if the 
involved relays independently choose which combinations of the 
messages to decode, there is no guarantee that the overall system 
of linear equations is solvable at the destination. In this article 
it is shown that, for a Gaussian fading channel model with two 
transmitters and two relays, always choosing the combination 
that maximizes the computation rate often leads to a case where 
the original messages cannot be recovered. It is further shown 
that by limiting the relays to select from carefully designed 
sets of equations, a solvable system can be guaranteed while 
maintaining high computation rates. The proposed method has a 
constant computational complexity and requires no information 
exchange between the relays. 

I. Introduction 

In wireless multiuser relay networks, both interference 
from multiple transmitters and noise degrade the system 
performance. To combat these issues, Nazer and Gastpar 
recently introduced a new relaying strategy called Compute- 
and-Forward (CaF) IT]. Their key idea is to decode an integer 
linear combination of the transmitted messages at intermediate 
relays, and then forward the combinations to the destination. 

Finding integer combinations that yield high transmission 
rates turns out, however, to be a complicated task. Particularly, 
finding the equation coefficients of the linear combinations that 
maximize the data transmission rate coincides with a Shortest 
Vector Problem (SVP) 121, for which various algorithms have 
been proposed a, a, 0. Unfortunately, these algorithms 
tend to be either highly complex or suboptimal. Recently, 
algorithms of polynomial complexity for finding the equation 
coefficients that maximize the rate have been presented fb], 
Q. However, choosing the coefficient vectors that maximize 
the instantaneous computation rate might result in an overall 
unsolvable system of linear equations at the destination due 
to linear dependency of the coefficient vectors - that is, the 
original messages might not necessarily be recoverable even 
if the combinations are successfully decoded at the relays. 

The message recoverability problem has been addressed in 
0, where precoding at the transmitters is used to increase 
the probability of receiving independent combinations at the 
destination, while Q, 11, El allow cooperation between 
the relays. These methods either require preprocessing at the 
transmitters or signaling between the relays. 

In Ea, this problem is mitigated by choosing a subset of 


relays with suitable equation coefficients. This approach can 
however not be used in the symmetric case, that is if the 
number of relays equals the number of transmitters. 

In this paper, we introduce a new, efficient approach to 
finding the coefficient vectors for a system with two trans¬ 
mitters and two relays. We observe that, in a Gaussian fading 
channel, there are only few coefficient vectors that typically 
maximize the computation rate. Based on this observation, 
we compile small coefficient vector candidate sets for both 
relays so that the probability of the relays choosing linearly 
dependent vectors vanishes. 

Our proposed method does not require jointly finding the 
coefficients at the relays, thus there is no need for cooperation. 
Furthermore, searching for appropriate coefficient vectors only 
over a small set of vectors reduces the computational complex¬ 
ity as opposed to solving the corresponding SVR 

The main benefits of our proposed method are twofold. 
Firstly, the end-to-end information outage probability vanishes 
as the message rate at the transmitters approaches zero. Sec¬ 
ondly, the computational complexity of our scheme is constant 
while still providing a relatively high throughput. Our findings 
are supported by extensive computer simulations. 

The paper is organized as follows. We give a brief introduc¬ 
tion of the CaF protocol in Section HI] and introduce methods 
for finding suitable equation coefficients in Section |ni| Sec¬ 
tion |IV| presents the performance metrics and corresponding 
numerical results, while Section |V] concludes the paper. 

II. The Compute-and-Forward Protocol 

In this article, we focus on a wireless multiple-access 
system with L = 2 transmitters and M = 2 relays, as 
illustrated in Figure [T] The first hop from the transmitters 
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Fig. 1. System model with two transmitters and two relays connected to 
a destination. The first hop is modeled as a wireless multiple access fading 
channel. The relays are connected to the destination with error-free bit pipes. 

to the relays is modeled as a Gaussian fading channel. The 
relays are connected to a destination with error-free bit pipes 










with unlimited capacities. The goal of the system is to reliably 
transfer information from both transmitters, via the relays, all 
the way to the destination. The relays apply the original CaF 
strategy introduced in IH, briefly exposed in the following. 

The transmitters want to communicate messages ' 

I = 1,2, where p is prime. Before transmiss 
messages are encoded into n-dimensional codewoi 
are subject to the power constraint ||x/|p < 
that throughout this paper, we assume that both tr 
use identical transmission powers, and that the re 
identical noise levels with normalized variance, i.e., I 
This is justified as allowing asymmetric SNR vali 
the performance metric used for simulations for 
considered strategies equally, and hence does not 
comparison results. 

Definition 1. The message rate at transmitter I is i 

nf = - loga p. 
n 

For the received signal at relay to, we use the following 
channel model: 

L 

ym — ^ ^ “f '^m^ (2) 

1=1 

where is additive white Gaussian noise with normalized 
variance — 1, and the channel coefficients are assumed 
to be i.i.d. and real-valued, hmi ^ ■A/’(0,1) IT]. Let h,„ = 
[hm .1 hm 2 ■■■ hmL deuote the channel vector for relay to. 

Channel state information is only available at the relays. 
Further, it is important to note that we assume that each relay 
only knows the channels to itself, i.e., relay 1 knows hu and 
hi 2 , while relay 2 knows ft, 2 i and ft, 22 - 

The key feature of CaF is to decode and forward linear com¬ 
binations of the transmitted messages. In this paper, we only 
consider integer combinations of the transmitted codewords. 

Remark 1. Ai noted above, each transmitter encodes its 
message Wi € into a lattice vector X;. Hence, although 
the destination ultimately wants to decode the messages w;, 
the relays decode linear equations involving the codewords 
X;, and it is thus meaningful to consider Ij-linear rather than 
¥p-linear combinations. The destination, upon reception of 
enough Z-linearly independent equations, can solve for the 
codewords and recover the original messages 

The integer combination of relay to is represented by an 
equation coefficient vector a^, more explicitely 

a^TT, ^ tlml flm2 ■" O’mL j . (3) 

These vectors form the equation coefficient matrix 

A = [ai a2 ■■■ aM f , (4) 

and the destination can recover the single codewords if and 
only if A is invertible, that is det(A) 0. 

Figure |2] shows that, if the relays always choose the equation 
coefficient vector that maximizes their instantaneous compu¬ 
tation rate, the probability of the matrix A being singular 


is noticeable. This phenomenon is more pronounced for low 
SNR values, and even more considerable for more than two 
relays and transmitters. In this article, we focus on the case 
M = L = 2. 



Fig. 2. Probability of the equation coefficient matrix A being singular for 
various numbers of relays M and transmitters L when independently choosing 
the coefficient vectors that maximize the computation rate. 


The most important performance metric of CaF is the 
so-called computation rate, whose meaning is explained in 
the following. If relay to achieves computation rate TZ^ = 
a^), then it is able to decode a linear combination 
of the codewords whose corresponding message rates satisfy 
< That is, relay to can support message rate 72.^ if 

7e?<min^, (5) 

where 7^® is the message rate at transmitter I, and Omi f- 0 
means that the message of transmitter I is included in the linear 
combination decoded by relay to. 

In the scenario considered in this paper, we have two 
transmitters with identical message rates = TZf = TZ^. 

Both relays can support this message rate if 

7e" < min{7^^,7^^}. (6) 


The main observations about the computation rate derived 
in H] are summarized in the following theorem. 

Theorem 1. In the above setup, a computation rate region of 


7e;,(h„,a„) = maxi log 
Z 


p 


+ P||ah„ - a™IP 


(7) 


is achievable. This expression is further maximized by choos¬ 
ing a to be 


OiMMSE — 


pyffi a 


1 + P||h„||2’ 

resulting in a computation rate region of 


77^(hm,a^) = - log”* 


P\K 


1-f P||h^||2 


( 8 ) 


(9) 


The computation rate is thus a function of the channel 
coefficients and the equation coefficients. While the channel 










coefficients are random variables, the relays are free to choose 
their desired equation coefficients. In this paper, we mainly 
focus on how to efficiently find equation coefficient vectors 
which yield high computation rates while ensuring that the 
vectors chosen by the relays are linearly independent. 


Proposition 1. The rotation U = [-?io] acts in a non¬ 
degenerate fashion on £. Any point [a] S £ has a unique 
pair [b] f [a] G £, so that U [a] = [b] and U [b] = [a]. 
Thus £ can be divided into disjoint sets £i and £2 so that 
£iU £2 = £, £ir ]£2 = 0, and \]£i = £^ 2 - 


III. Selecting the Equation Coefeicient Vectors 

In this section we present three methods in detail for finding 
equation coefficient vectors a^ at the relays, all of which 
require no cooperation between the relays. 

The main interest in this article is to guarantee a solvable 
system of equations at the destination, while still being able 
to support high message rates. To that end, we will include 
the following indicator function in our performance metric: 


J 1 if det(A) f 0 , 
{det(A)/o} - I Q det(A) = 0. 


( 10 ) 


1. Exhaustive search: Each relay searches over all possible 
equation coefficient vectors and chooses the one that maxi¬ 
mizes its instantaneous computation rate (|9]l. 


Lemma 1. ^ Finding the coefficient vector that maximizes 
the instantaneous computation rate is equivalent to solving 


a.m = argmin a^Ga, (11) 

aeZJ:'\{(0,...,0)} 


where 


G = Il- 


1 + P||h„|| 2 ’ 


( 12 ) 


and corresponds to finding the shortest vector in the lattice 
whose Gram matrix is G. 


While maximizing the computation rate individually at each 
relay yields the highest message rate supported by the given 
channel realizations, there is no guarantee that the chosen 
vectors form an invertible matrix at the destination. 

2. Proposed splitting method: In our proposed method, we 
split possible coefficient vectors into two disjoint sets Vi and 
V 2 . Set Vi is pre-assigned to one relay, while set V 2 is pre¬ 
assigned to the other. Each relay independently chooses the 
vector within its set that maximizes the computation rate (|9l). 

Remark 2. Let am = [“>ni “^ 2 ]^ be a solution to (fTTI) . Note 
that the matrix G is symmetric, that is results in 

T T [ X z •[ 2 .0 i2 

[ 2 : y J 2(lml0^m2Z H“ ^rn2y' 

Thus if [“ml 0 ^ 2 ]^ is a solution to (|9]), so is [ 0 'm 2 j 

In particular, we can fix ami > 0. 

One may thus divide the set of candidate vectors 
Z^\{(0, 0)} into equivalence classes modulo the Z 2 -acf;on 
of ±1. Solutions to (HB are sought for in the quotient 
space £ = (Z^\{(0, 0)}) mod Z 2 . The elements in £ are 
equivalence classes [a] of vectors a G Z^\{(0,0)} up to 
Z 2 , i.e. pairs of vectors, [a] = {a, —a}. A unique way of 
representing elements in £ is by points in Z^\{(0,0)} in the 
right half plane, including the upper y-axis, but not the lower 
part, i.e., with vectors (ai, 02 ), oi > 0 or ai = 0, 02 > 0. 


Proof: Consider the vector a = (01,02) G Z^\{(0, 0)}. 
Its equivalence class is {(oi, 02), (—01, —02)}. The action of 
U on a vector commutes with scalar multiplication by ±1. 
Thus U[a] = {(02,01), (-02,01)} = [b] G £. Clearly, [a] 7 ^ 
[b], as a 7^ (0, 0), and U [b] = {(-oi, -02), (01,02)} = [a]. 
The action of U thus divides £ into a disjoint set of pairs. 
One element of each pair can be taken to £ 1 , the other to £ 2 - 


Eollowing Remark m we define the following set. 

V = {{x,y) €l?\x>l,y + y^ < 

d G Z\{0},gcd(x,y) = 1} . 

Let Vi,V 2 C V such that |Vi| = IV 2 I, Vi U V 2 = V and 
El n = 0 , and define 

1/1 = ViU 1(0,1)}, E 2 = V 2 U 1(1,0)}. (14) 


While preventing the relays from selecting arbitrary coeffi¬ 
cient vectors decreases the expected computation rate at the 
relays, choosing linearly independent vectors ensures that the 
messages can be recovered at the destination - provided that 
the relays are able to support the message rate at the transmit¬ 
ters. This claim is justified in the following proposition. 

Proposition 2. For any choice of sets Vi, V 2 as proposed 
above, we have almost surely 


min{7^},7^B • l{det(A)#o} > 0- (15) 

Proof: Since Vi D V2 = 0 , we have l{det(A)/o} = 1 - 
It hence suffices to show that there exist vectors ai G Vi, 
a2 G V2 such that for every channel realization h = [ tii /12 ]^, 


2 


P\h^am\^ Y\ 

i+piihip; 


or equivalently 0 < ||am ||2 - < 1 - 

Choose ai = [10]^ G Vi, a 2 = [01]^ G V 2 . Then, for 
m = 1,2 we have ||am|P = 1 , and, since 


Ih^£ 


hf if TO = 2 , 
if TO = 1, 


for TO = 1,2 we have 0 < 
0 < 1 - 


Phi 


Phi 


< 1 , hence 


1 + Pihl + hl) - 


as required, where equality holds if and only if hm = 0 - ■ 

The following proposition gives a further criterion for 
designing the sets of vectors Vi and V 2 - 


Proposition 3. Let Vi and V 2 be two sets of vectors, chosen 
as above, and denote by the expected computation 







rate of the first and second relay using the assigned sets, 
respectively. If there exists a rotation matrix U such that 

V:Vi^V 2 ; at-^Ua (16) 


is a bijective isometry, we have 


T^v,=T^v.- 


(17) 


Proof: Assume U is such a matrix. Then, for every a G 
V 2 we find a G Vi such that a = Ua. Moreover, ||a|p = 

liu&ip. 

Let a = Ua, h = U^h. Then, for a fixed power P, and 
since U^U = I 2 , 


7e’'(h,a) = ilog+ 

= ilog- 



1 + P||h|p; ) 



p|h'^ua|2\ 
T+p\\W) ) 


= iiog" 
= iiog- 



P|(Uh)^Uap \ 
l + P||Uh||2j j 



Pjh^ap \ \ 

^ + P\\w) ) 


= P’'(h,a). 


It follows P’'(h,a) = 7e’'(Uh,Ua) = P'-(h,a). 

To conclude the proof, note that the distribution of t 
channel h is rotation invariant. We thus have 


= Uh 

max7^’'(h, a) 

aeVi 

= Uh 

max 7eUh,Ua) 
u&eVi 

= Uh 

max7^’'(Uh,Ua) 

_aey2 

= 

Pi. 

max7^'’(h, a) 
aeVG 

II 







3. Multiple Access: For the sake of comparison, we also 
consider a simple multiple access method where each relay 
decodes the message corresponding to the stronger channel 
coefficient while treating the other message as noise. More 
concretely, if = [hmi hm 2 ]^ is the channel observed by 
relay m, then the relay decodes Xi iff hmi > hm 2 , and 
decodes X 2 otherwise. In this scenario, for i,j = 1,2, i j, 
if \hmi\ > \hmj\, message x^, can be decoded by relay m if 


K < 2 log 



P\hrm\^ \ 

1 + P\hmj P / 


(18) 


We say that a determinant error occurs if both relays decode 
and forward the same message. 


IV. Simulation Results 

In this section, we present extensive simulation results for 
the end-to-end performance of the methods introduced in 
Section [111] for the 2-by-2 relaying system described earlier. 

We start by defining two sets of vectors Vi and V 2 as in 
(O for the proposed strategy. Note that for a given channel 


realization hm, the equation coefficient vector maximizing 
the computation rate at relay m is the one that best aligns with 
hm IB. It is thus crucial to provide each relay with vectors 
that can approximate channel vectors lying in any direction. 

Further, to achieve low computation complexity, we con¬ 
struct disjoint sets of small cardinality. For the following 
simulations, we fix d = 4 in Cl, resulting in sets of 
cardinality \V\ = 14, \Vi \ = IV 2 I = 8. 

Proposition 4. Splitting as discussed above exists, where V 2 = 
UVi. For such splittings, min{Pj^,P2} ’ l{det(A)} > 0 for 
any channel realization, almost surely. Further, TZy^ = Pyy 

Proof: Choose Vi such that Va G Vi, either [a] G or 
U [a] G £ 1 , but not both. Then according to Prop. [T] UVi G 
£ 2 , and, Vi n V 2 = 0- Note that if [(1,0)] G £ 1 , we have 
U [(1, 0)] = [(0,1)] G £ 2 - Thus, according to Prop.|2] the first 
statement holds and by Prop. [3 the second statement follows. 

■ 

One example of such Vi and V 2 is depicted in Figure [3] 



Fig. 3. All vectors that maximize the computation rate for 10® channel 
realizations (left) and sets of vectors Vi, V 2 for the proposed method (right). 


We now compare the three introduced methods. Our two 
main performance metrics are minimum end-to-end rate and 
end-to-end outage as described in the following. 

Minimum end-to-end rate: This metric is defined to be the 
expected value of the minimum of the two computation rates 
at the relays. If the coefficient matrix is not invertible, this 
metric is defined to be zero as the destination cannot recover 
both of the original messages. The minimum end-to-end rate 
can be expressed as 

nQ = E (min {TZ[, 7^^} • l{det(A)#o}) , (19) 

where E {■) denotes the expected value, and is the 
computation rate at relay m defined as in (|9]l. We illustrate 
this performance metric for varying SNR values in Figure H 
The proposed splitting method performs relatively close to 
the exhaustive search method, and clearly beats the multiple 
access strategy (cf. Section 11111-31. It is to be expected that 
the splitting method cannot outperform the exhaustive search 
whenever it has positive rate (cf. (fT^ l in terms of end-to-end 
rate, since the exhaustive search has a considerably higher 
number of coefficient vectors from which to choose. The 























Fig. 4. Expected minimum end-to-end rate vs. SNR. 


drawback of the exhaustive search is, however, its high com¬ 
plexity compared to the constant complexity of our proposed 
method. We argue that, in many cases, it is more important 
to have a fast, lightweight algorithm that finds relatively high 
rates, rather than strictly maximizing the rate at the cost of 
computation time and power. 

End-to-end outage: The system is said to be in outage if at 
least one of the relays cannot support the symmetric message 
rate at the transmitters, or if the equation coefficient matrix is 
not invertible. The outage probability with symmetric message 
rates TZ^ can be expressed as 

p = Pr((min{7^^,7^^} < 7e^) U (det(A) = 0)). (20) 

Figure 13 presents simulated end-to-end outage probabilities 
for a fixed SNR value as a function of the symmetric message 
rates. While for the exhaustive and multiple-access strategies 



Fig. 5. End-to-end outage as a function of the message rate for SNR 10 dB. 

the outage probability converges to a non-zero value for 
decreasing message rate, for the proposed splitting method 
it goes to zero as the message rate goes to zero, as indicated 
by Prop. 121 

V. Conclusions and Future Work 
I n this article it was shown that pre-assigning sets of 
equation coefficient vectors to intermediate relays in the 


framework of Compute-and-Forward relaying helps preventing 
the overall system of linear equations at the destination from 
being singular, while maintaining high computation rates and 
thus being able to support high message rates. The advantage 
of the proposed method its constant complexity, as opposed 
to solving a hard shortest vector problem. Moreover, no 
cooperation between the relays is required. 

In this article, we restricted ourselves to real-valued chan¬ 
nels, as well as the case of only two transmitters and relays. 
It was however observed that the problem of the coefficient 
matrix being singular is even more dramatic for larger number 
of transmitters and relays. Therefore, as a natural extension, 
future work includes generalizing the introduced method for 
an arbitrary, not necessarily symmetric number of transmitters 
and relays, as well as considering complex-valued channels. 

Moreover, while design criteria for the sets of vectors have 
been mentioned, it still remains open to describe analytically 
a way of finding the best possible disjoint partition of the 
equation coefficient vectors. 
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